ENTROPY PRODUCTION IN QUANTUM SPIN SYSTEMS. 



by David Ruelle*. 



Abstract. We consider a quantum spin system consisting of a 
[ finite subsystem connected to infinite reservoirs at different tem- 

peratures. In this setup we define nonequilibrium steady states 
■ and prove that the rate of entropy production in such states is 

. nonnegative. 
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For several decades, Joel Lebowitz has been the soul of research in statistical me- 
chanics. He now plays a central role in the development of new ideas which reshape 
our understanding of nonequilibrium. The present paper, dedicated to Joel on his 70-th 
birthday, extends some of the new ideas to quantum systems. 

Introduction. 

Consider a physical situation where a "small" system S is connected to different 
"large" heat reservoirs R a (a = 1,2,...) at different inverse temperatures (3 a . We want 
to define nonequilibrium steady states for the total system L = S + Ri + R2 + . . ., and 
verify that the rate of entropy production in such states is > 0. The model which we 
discuss in this paper is that of a fairly realistic quantum spin system. In what follows 
we first describe the model and state our assumptions (Al), (A2), (A3). In this setup we 
introduce nonequilibrium steady states p as states which, in the distant past, described 
noninteracting reservoirs at different temperatures. Under suitable conditions we check 
that our definition does not depend on where we place the boundary between the small 
system and the reservoirs. Our definition of the entropy production e p also does not 
depend on where the boundary between the small system and the reservoirs is placed. 
With this definition we prove e p > 0. By contrast with an earlier paper [4], we omit here 
assumptions of asymptotic abelianness in time which are difficult to verify, the definition 
of nonequilibrium steady states is more general, but we obtain less specific results. 

Description of the model.* 

Let L be a countably infinite set. For each x G L, let 7i x be a finite dimensional 
complex Hilbert space, and write Hx = ®iex^i if X is a finite subset of L. We let Ax 
be the C*-algebra of bounded operators on Hx, and if Y C X we identify Ay with a 
subalgebra of Ax by the map Ay Ay <S> 1h x \y c Ax- We write L as a finite union 
L = U a >oR a , where Rq = S is finite (small system) and the R a with a > are infinite 
(reservoirs) . We can then define the quasilocal C* algebras A a , A as the norm closures of 

U Ax , |J A x 

XcR a XCL 

repectively. Note that all these algebras have a common unit element 1. In this setup we 
assume that an interaction $:Ih ®(X) is given such that $(X) is a selfadjoint element 
of Ax for every finite X C L. Also, for each reservoir, we prescribe an inverse temperature 
Pa > and a state a a on A a - 

The assumptions (Al), (A2), (A3). 

(Al) The interaction $ satishes 

II$IIa = $>™ a su P E ||*p0||<oo 

n>0 xEL X3x:cardX=n+l 

* See [3], [1]. 
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for some A > 0. 

The importance of this assumption is that it allows us to equip A with a one-parameter 
group (a*) of automorphisms* defining a time evolution. Introduce a linear operator 5 : 
UxclAx — ► A such that 

6A = i J2 if AeAx 

Y:YnX^<D 

If A E Ax, one checks that 

\\5 m A\\ < ||A||e AcardX m!(2A- 1 ||$|| A ) m 
The strongly continuous one-parameter group (a*) of *-automorphisms of A is given by 

a 1 A = Y —d m A 
' m! 

m=0 

if A G Uxc-L-4x and \t\ < A/2||$||>. (More generally one could take A e where ^ is 
defined in the Appendix). Let 

XcA 

for finite Ac I and A G i. Writing A — > L if A eventually contains each finite A C L, 
we have then 

lim \\e ttHA Ae- UHA - a 1 A\ \ = 

A^L 

uniformly for t in compact intervals of R. 

(A2) $(A) = if A n S = 0, A n R a + 0, A n R b ^ for different a,b>0. 

Note that the description of the interaction $ is somewhat ambiguous because any- 
thing ascribed to $(A) might also be ascribed to <&(Y) for Y D X . Condition (A2) means 
that in our accounting, if a part of the interaction connects two different reservoirs, it must 
also involve the small system S. 

(A3) If a > 0, let $ a be the restriction of the interaction $ to subsets of R a and write 

H aA = = ff fl«nA 

xcR a r\A 

Let also the interactions ^(a) ° e given such that 

||*(A)|U< K<oc (1) 

* See [1] Theorem 6.2.4 (or [3] Section 7.6). 
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and write 

B aA = Yl *(A)W 

xci? a nA 

We assume that, for a suitable sequence A — >■ L, 

Tr-u ( e -Pa(H aA +B aA )A\ 
hm — — — 2 -3-777 — -= — r — = a a (A) 

if A G A a : this dehnes a state cr a on A a , depending on the choice of (^(a)) and the 
sequence A — > L. Furthermore we assume that for each Gnite X there is Ax such that 
*(A)C*0 = if A D A x and 7cX; therefore 

||[BaA,4||=0 (2) 

if A D Ax and A G .Ax- 
In particular we can take all ^(a) = 0. Using (3) below, it is readily verified that 
a a is a /3 a -KMS state (see [2]) for the one-parameter group (a*) of automorphisms of «4 a 
corresponding to the interaction $ a . [I do not know which of the /3 a -KMS states can be 
obtained in this manner]. 

Note that the assumptions (Al), (A2), (A3) can be explicitly verified in specific cases. 
From (A3) we obtain the following result. 

Lemma. 



lim We^iHaA+B^) Ae -zt{H aA +B aA ) _ & t A u = q / 3 x 

A — >L a 

for a > 0, and 

lim | \ e U{HA+ ^ a >° B ^Ae- U{HA+ ^«>o b « a) -a t A\\ = (4) 

A^L 

uniformly for t in compact intervals of R. 

We prove (4). Write a A A = e lt{HA+ ^«>o BaA) ie" !t(ffA+ ^«>» BaA) and 5 A A = i[H A + 
E«>o S «A' A \- If A G U X A X we see using (1) that 



m=0 



converges uniformly in A for \t\ < A/2(||<E>|| A + K). Using also (2), it is shown in the 
Appendix that 5™ A — > 5 m A in A when A — > L. Therefore 

lim \\a\A - ofAW = 

A— >L 



when A e UxAx, uniformly for \t\ < T < A/2(||$|| A + K). But the condition A e U X A X 
is removed by density, and the condition \t\ < T < \/2(\\<&\\\ + K) by use of the group 
property. The proof of (3) is similar, [j 

The KMS state a. 

The interaction J2 a >o Pa® a, evaluated at X is f3 a <& a (X) if X C R a and if X is not 
contained in one of the R a - The corresponding one-parameter group (/?*) of automorphisms 
of A has, according to (A3), the KMS state* a = ® a >o^a where a"o is the normalized trace 
on Ao = As- In fact 

„( a\ _ ]]m Tr HA (eM-j:aPa(HaA + B aA ))A) 



Nonequilibrium steady states. 

We call nonequilibrium steady states (NESS) associated with a the limits when T — > oo 

of 

tJ dt{atya 

using the w*-topology on the dual A* of A. With respect to this topology, the set E of 
NESS is compact, nonempty, and the elements of E are (a*)*-invariant states on A. 

This definition generalizes that given in [4] where, under stringent asymptotic abelian- 
nes conditions, the existence of a single NESS was obtained. 

Dependence on the decomposition L = S + R\ + R2 + . . .** 

Our definition of a, and therefore of E depends on the choice of a decomposition of 
L into small system and reservoirs. If S is replaced by a finite set S' D S and the R a by 
correspondingly smaller sets R' a C R a one checks that (Al), (A2), (A3) remain valid. If 
$^ is the restriction of $ to subsets of R' a , the replacement of f3 a & a by ^ f3 a &' a changes 
(/?*) to a one-parameter group (/?'*) and cr to a state cr'. These changes are in fact bounded 
perturbations covered by Theorem 5.4.4 and Corollary 5.4.5 of [1]. The map a — > a' (of 
KMS states for (/?*) to KMS states for (/?'*)) is nonlinear (as can be guessed from (5)) and 
therefore we cannot expect that ^ dt(a t )*a' has the same limit as ^ J Q T (it (a*)*cr in 
general, but the deviation is not really bad. The (central) decomposition of KMS states 
into extremal KMS states gives factor states. If a is assumed to be a factor state, and 
(a*) is asymptotically abelian, one finds that lim ^ J Q dt(a t )*a does not depend on the 
decomposition L = S + Ri + R2 + ■ ■ ., as the following result indicates. 

Proposition. 

* The state a corresponds to the inverse temperature +1 rather than the inverse tem- 
perature — 1 favored in the mathematical literature. 

** This section and the following Proposition are in the nature of a technical digression, 
and may be omitted by the reader essentially interested in the positivity of the entropy 
production. 
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Using the above notation, assume that a is a factor state, and that 

lim ||[a*A,B]|| = 

when A, B E A. Then, when T — > oo, 



lim^ I eft (a*) V = lim ^ / dt(a z )*a 
T Jo T J 



Let us introduce the GNS representation (H, 7r, O) associated with cr so that if 

1 r T 

p = lim— / .it (a*)*cr 

-t ./n 



we have 

p(A) = lim 1 / dt(n,7r(a*^)n) 
io 

By restricting T to a subsequence we may assume that in the weak operator topology 

lim ^ f dt ilia 1 A) = A E n{A)" 
T Jo 

and by assumption we also have A e 7r(„4)', hence A e 7r(^4)' n 7r(.A)" = {Al} since a is a 
factor state. 

But we may write cr'(-) = (O', 7r(-)0'): this follows from the perturbation theory of [1] 
(see proof of Theorem 5.4.4). We have thus 

lim^ / dt a' {a* A) = lim ^ I dt (fi' >(a* A)fl') 
T Jo T Jo 

= lim^ I dt (0, 7r(a*A)n) = lim 1 / cfta(a*A) 

as announced, [j 

Entropy production. 

For finite A C L we have defined 

XcA 

but Hl-, Hn a do not make sense. We can however define 

[H L ,H Ra ] = hm[H A ,H RanA ] = lim[H A ,H aA ] 

A^L A^L 



We have indeed 

[Ha, H aA ] = [H A - H aA , H aA ] = [H A - ^ H 6A , H aA ] 

b>0 

and (A2) gives 

b>0 xES X:xeXcA v ' 

[implying the existence of the limit limA->z,(HA - J2b>o H b\) = H L - J2b>o H Rb e ^]. 
Using (Al) we obtain 

\\[<f>(X),H aA }\\ < 2A- 1 ||$|| A ||$(X)||e AcardX 

hence 

\\MX),H aA ]\\ <2A- 1 ||$|| A e A ||$|| A 

XBx 

and [Ha, H aA ] has a limit [H L , H Ra ] e A when A — > L with 

||[H L ,H i? J||<2card^A- 1 e A ||$||| 

The operator 

*[#L,H fl J 

may be interpreted as the rate of increase of the energy of the reservoir R a or (since 
this energy is infinite) rather the rate of transfer of energy to R a from the rest of the 
system. According to conventional wisdom we define the rate of entropy production in an 
(a*)*-mvariant state p as 

e p = J2Pa P (i[H L ,H Ra }) 

a>0 

(this definition does not require that p G E). 
Remark. 

If we replace S by a finite set 5" D S and the R a by the correspondingly smaller sets 
R' a C -R a , we have noted earlier that (Al), (A2), (A3) remain satisfied. As a consequence 
of (Al) we have 

i[H L , H Ra - H K ] = lim i[H A , H aA - H' aA ] = lim ${H aA - H' aA ) 

A. — ^/j TV — ^i/ 

(where the operator 5 has been defined just after (A3)), hence 

p(i[H L ,H Ra -H K })= hm L p(5(H aA -H' aA ))=0 

i.e., the rate of entropy production is unchanged when S and the R a are replaced by S' 
and the R' a . The reason why we do not have p{i[Hl, HrJ) = is mathematically because 



H Ra is "infinite" (H Ra ^ A), and physically because our definition of p(i[H L , H Ra ]) takes 
into account the flux of energy into R a from S, but not the flux at infinity. 

Theorem. 

The entropy production in a NESS is nonnegative, i.e., e p > if p E E. 
We have seen that 

[H L ,H Ra ]= lim [H A ,H aA ] 

A^L 



= lim [H A - V" H bA , H aA ] 

A^L L — * 

Therefore, using (A3) and [H bA + B bA , J2 a>0 (H aA + B aA )} = 0, we find 



ft a [Hl 7 H Ra ] = |m[ffA - £ # b A, £ A,#oa] 
a>0 " b>0 a>0 

= lim [if A - V H bA , V /3 a (tf aA + S oA )] 
a — > r, g * g * 



A — >.L 

&>0 a>0 

= Um [# A + Bh ^ ^(HaA + BaA)] 
6>0 a>0 

in the sense of norm convergence. 

We also have, for some sequence of values of T tending to infinity and all A E A, 

p(A) = lim 1 / dtaia'A) = lim lim 1 / rfta(a A A) 

where, by (4), 

a A A = e v AT ^a> oA 'Ae v A ^z^a>o aA ' _ > a z A m norm 

when A — > L, uniformly for f G [0, T]. 
Write 

i?BA = #A + ^ B aA 

a>0 

G A = ^/3 a (tf aA + S aA ) + logTr HA exp(- ^/3 a (tf aA + B aA )) 

a>0 a>0 

Then the entropy production is 

e p = p{iY j (3 a [H L ,H Ra ])= lim lim 1 / df a(e^ A [ff BAj G A ]e"^ A ) 



and the convergence when A — > L of the operator (e ltHBA [LTba, G r A]e _ltHsA ) is uniform for 
t E [0, T]. According to (A3) we may choose the A tending to L such that Tr nA e~ GA (•) 
tends to cr(-) in the w*-topology, hence 

e p = lim lim-/ Tr WA (e- GA e^ SA [if BA , G A ]e- itHBA ) 

T— »oo A— »L _Z 

= lim lim 1 / dtTr nA (e- G ^(e UH ^G A e- UH ^)) 

T— »oo A^L 1 Jq OX 

= lim lim i(Tr WA (e- GA e iTHflA G A e- iTHBA )-Tr WA (e- GA GA)) 

T — >oo A — >L 1 

and the Theorem follows from the Lemma below, applied with A = Ga, U = e tTHBA and 
0( s ) = -e" s . 

Lemma.* 

Let A, U be a hermitean and a unitary n x n matrix respectively, and ^ : R ^ R be 
an increasing function. Then 

tr{(j){A)U AU- 1 ) < tr((j)(A)A) 
By adding a constant we may assume <j> > 0. We write the spectral decomposition of 

A as 

n n 

A = a J^' = aiPi + J2( a J - a J-^ p j 

3 = 1 J=2 

where ai < . . . < a n are the eigenvalue of A, Ej the spectral projections, and Pj = 
Ej + . . . + E n . Therefore 

n 

0(A) = 0(a 1 )P 1 + - ttaj-^Pj 

i=2 

and 

tr(0(A)^^ _1 ) = 0(ai)tr(P 1 t/^C/- 1 ) + ^(0(a^) - 0(a J -_i))tr(P J -^^- 1 ) 

J=2 



We have 



* As R. Seiler kindly pointed out to me, this lemma can be obtained readily from O. 
Klein's inequality 

tr(f(B)-f(A)-(B-A)f'(A))>0 
where A, B are hermitean and / convex: take B = U AXJ~ X and <j) = f. 
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where c kl = tr (E k UEiU x ) > and J2k=i c ki = h Ya=i c ki = 1 - Hence 

n n n 
1=1 k=j l=j 

because < Yl=j c kl < 1> Ya=i Yl=j c ki = n - j + 1 and a-y < . . . < a n . Finally 

n n n n 

tiMAlUAU- 1 ) < 0(ai) a i + £(0( a j) - Haj-i)) J2 ai = Yl ^ a j) a J = *■*(<!>( A) A) 

1=1 j=2 l=j 3-1 

proving the Lemma, [j 
Remark. 

We have 

a>0 

because 

-Y,P(i[H L ,H Ra }) = ]hnp(i[H A ,H A - ^^a]) = ^p(a* £ *(*))|t=o = 

a>0 ' a>0 "' X:XnS#0 

where we have used the fact that p is (a*)* -invariant. In particular, in the case of two 
reservoirs 

0<e p =(p 1 -p 2 )p{i[H L ,H Rl }) 

so that if the temperature P^ 1 is less than Z^ 1 , i.e., fii — f3 2 > 0, the flux of energy into 
Ri is > 0: heat flows from the hot reservoir to the cold reservoir. 

Proving strict positivity of e p . 

It is an obvious challenge to prove that e p ^ 0. A natural situation to discuss would 
correspond to R a = Z" and $ a translationally invariant. But we need then v > 3 as 
discussed in [4]. Indeed, for v < 3 one expects a nonequilibrium steady state to be in fact 
an equilibrium state at a temperature intermediate between the original temperatures of 
the reservoirs. Instead of a quantum spin system as described above, a gas of noninteracting 
fermions would probably be easier to treat first. 

Appendix: the algebras A\. 

The purpose of this Appendix is to complete the proof of (4) by establishing (10) 
below. 

If A > 0, a norm \\.\\\ is defined on UxAx by 

||A|| A = inf^||A x ||e AcardX 
x 
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where the inf is taken over finite decompositions A = J2x Ax with A x G Ax- We let A \ 
be the completion of UxAx with respect to that norm. If A > \i > there are natural 
continuous maps A\ — > A^ — > A, which are not claimed to be injective. 

If | |<E>| | \ < oo and A x G Ax the formula 

SA x = i M Y )' A *] 

Y:YnX^0 

defines an element of A\. If A > /i > 0, and \\&\\\ < oo, one also checks that 5 m defines a 
map A\ — > wA^j such that 

||^IU< ||A|| A m!(2(A-^)- 1 ||$|Ur (6) 
We have 5 A = 5 A + 5 A , where 

5 A A = z[tf A ,A] , ^A = z[^S aA ,A] 

a>0 

and (1) and (6) (for m = 1) yield 

||<L4|| M <||,4|| A .2(A-/x)- 1 p|| A 

||<MU< II^IU.2(A-^)- 1 ||$|U 

||^IU< ||-4_!U-2(A — a*)- 1 ^ 

Given e > and A e A\ we can find X such that A = A\ + A 2 with Ai G and 
||A 2 ||a < €. Therefore 

-6 A )A\\v< \\(S-S A )A 1 \\ li + \\SA 2 \\ li + \\6 , A A 2 \\ li + KA 2 \\ li 

- <y A )^ilU + e.2(A - ^)- 1 (2||$|| A + K) (7) 



Taking A D Ax we also have 

5^i = 

by (2), and 

(5-S' A )A 1 =i Yl Mn^i] 

Y:Y<£A,YnX=£® 

so that 

||(5 - 5 A )^i|U < II^i||a.2(A - ^^II^IIxa (8) 

where ||$'|Ua = sup x6X Y,ybx y<£x e {caTdY ~ 1)X \\$(Y)\\. When A -> L we have ||$'|Ua -»• 
and (7), (8) yield 

Urn ||(*-* A )4| M = (9) 
11 



We can now prove that, if ||$||a < °o and A e A\, 



lim \\8 m A-6%A\\ = (10) 



We have indeed 

and, using (6), 

hence, by (9), 
so that, using (6), 



m— 1 



k=0 

\s k A\\ 2X/3 <\\A\\ X M(j\m x y 

lim \\(8-8 A )8 k A\\ x/3 = 

A. — * L 



\\5 A a - k - 1 (5-5 A )5 k A\\ < \\8^- k - 1 (8-8 A )8 k A\\o 

< - 5 A )5 fc A|| A/3 (m - k - l)!(^p|| A r- fc - 1 

which tends to zero when A — > L. This concludes the proof of (10). 
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